Thermodynamic properties of a dipolar Fermi gas 
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Based on the semiclassical theory, we investigate the thermodynamic properties of a dipolar 
Fermi gas. Through a self-consistent procedure, we numerically obtain the phase-space distribution 
function at finite temperature. We show that the deformations in both momentum and real space 
become smaller and smaller as the temperature is increased. For the homogeneous case, we also 
calculate pressure, entropy, and heat capacity. In particular, at the low-temperature limit and in the 
weak interaction regime, we obtain an analytic expression for the entropy which agrees qualitatively 
with our numerical result. The stability of a trapped gas at finite temperature is also explored. 
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I. INTRODUCTION 
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The experimental success in creating ultracold 
^°K*^Rb molecular gas near quantum degeneracy [ll-0| 
has drawn considerable attention to in studying the fun- 
damental properties of degenerate dipolar Fermi gases. 
Within the framework of the semi-classical theory, the 
ground state properties, the collective excitation, and the 
free expansion dynamics of a normal state dipolar Fermi 
gas were studied theoretically [5 -7]. A recent theoret- 
ical work based on a variational approach reveals that, 
due to the Fock exchange interaction, the momentum 
distribution is stretched along the direction of dipole mo- 
ment such that the Fermi surface becomes an ellipsoid [sj . 
This result was confirmed numerically for both homoge- 
neous [9] and trapped systems. Taking into account 
the effect of the exchange interaction, further theoretical 
work regarding the normal state of the zero tempera- 
ture dipolar Fermi gas includes studying the free expan- 
sion [3, [ll| , collective excitation 
and the Fermi liquid properties [ij 

Another interesting feature of the dipole-dipole inter- 
action is that the partially attractive dipolar force is 
responsible for the formation of anisotropic BCS pair- 
ing fls'-'lG'l. With the recent experimental development 
on control the hyperfine states of ^"K^^Rb molecules [1, 
0], the BCS pairings in a mixture of fermionic polar 
molecules with two different hyperfine states are also 
studied theoretically (T7l - [l9| . In particular, the effects 
of the Fock exchange interaction to pairing was consid- 
ered in Ref. 

In this paper, we extend our previous work on the 
ground state properties of the dipolar Fermi gases to fi- 
nite temperature case. Employing the semi-classical the- 
ory, we numerically obtain the phase-space distribution 
function through a self-consistent procedure. We show 
that the deformations of the distribution function in both 
momentum and real space become smaller and smaller as 
the temperature is increased. For homogeneous system, 
we also calculate the thermodynamic quantities such as 
pressure, entropy, and heat capacity. In particular, at low 
temperature limit and in weak interaction regime, we de- 



rive an analytic expression for the entropy, which agrees 
qualitatively with our numerical result. For the trapped 
gases, we also explore the temperature dependence of the 
stability. 

The remainder of this paper is organized as follows. 
In Sec. |TT1 we introduce our model and briefiy outline 
the semi-classical theory for a dipolar Fermi gas at finite 
temperature. The numerical and analytical results are 
presented in Sec. IIIIl Finally, we conclude in Sec. |IVl 



II. THEORY 

Here we consider a system of N spin polarized dipolar 
fermions interacting via dipole-dipole interaction 
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where Cd — (P/{4:Treo) with d being the electric dipole 
moment. For simplicity, we have assumed that the 
dipole moments of all fermions are orientated along z- 
axis. Additionally, we shall consider both homogeneous 
and trapped systems. In the former case, C/cxt can be 
regarded as a box of volume V; while for the latter, the 
trap is assumed to be a harmonic potential with axial 
symmetry, i.e., 

C/e.t(r) = \mZj'\-'/'ix' +y' + X'z'), (2) 

where ZJ is the geometric average of the trap frequencies 
and A is trap aspect ratio. 

Within the framework of the semiclassical theory, the 
thermodynamic properties are completely characterized 
by the phase-space distribution function /(r,k), which, 
at finite temperature T, satisfies the Fermi-Dirac statis- 
tics 



/(r,k) 



g(£(r,k)-Ai)/fcBT 
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where fi is the chemical potential introduced to fix the 
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FIG. 1: (Color online), (a) and (b) represent, respectively, 
the phase-space distribution functions /(fcp, 0) and /(O, fcz) for 
Dji = 0.3. The corresponding temperatures are ksT/e'p — 1 
(solid lines), 0.5 (dashed lines), 1/3 (dash-dotted lines), and 
0.1 (dotted line), (c) The temperature dependence of a for 
various dipolar interaction strengths. The horizontal lines 
denote the values of a at zero temperature. 
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FIG. 2: (Color online). The temperature dependence of i5kin 
(a), -Ecxc (b), and \i (c). The dimensionless dipolar interaction 
strength are = (dash-dotted lines), 0.2 (dashed lines), 
and 0.3 (solid lines). 



total number of particles such that 

drdk , 

-/(r,k) = 7V. 
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The dispersion relation of the quasi-particle takes the 
form 
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f{r',k')Vd{r-r') 



dk' 



■/(r,k')Fd(k-k'), 



(5) 



where Vd(k) = — 1^ is the Fourier transform 

of Vd{r) and the last two terms represent the mean- field 
potentials originating from Hartree direct and Fock ex- 
change interactions, respectively. We remark that the lo- 
cal density approximation has been employed to obtain 
Eq. (|3]) for the trapped system. 

Equations ([3|)-(l5]) form a closed system of equations 
which can be solved numerically through an iterative pro- 
cedure to obtain the phase-space distribution function 
/(r,k). Within semi-classical approximation, the total 
energy at the finite temperature takes the same form as 
that at zero temperature, i.e.. 
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drdkdr'dk7(r, k)/(r', k') Vd(r - r' 
drdkdk7(r,k)/(r,k')Frf(k - k'), 



where the first line represents the kinetic (ii'kin) and po- 
tential {E-pot) energies, and last two lines are, respec- 
tively, the direct (i?dir) and exchange (i?cxc) interaction 
energies. Finally, we point out that, for axially sym- 
metric system, the phase-space distribution function re- 
duces to /(r, k) — f{p,z,kp,kz), where p — \/x^ + y'^ 

and kp — y^fc^ -I- ky . This fact can be used to simplify 
the numerical integration T^. 



III. RESULTS 

In this section, we present our results on the thermo- 
dynamic properties of a dipolar Fermi gas. For a homo- 
geneous gas, we first obtain the phase-space distribution 
function numerically, which allows us to calculate other 
thermodynamic quantities, such as pressure, entropy, and 
heat capacity. At low temperature limit and in weak in- 
teraction regime, we also derive an analytic expression 
for the entropy, which qualitatively agrees with our nu- 
merical result. In the second part of this section, we 
present the numerical results on the thermodynamics of 
the trapped systems. 



A. Homogeneous case 

For a homogeneous gas, the phase-space distribu- 
tion function reduces to a function of momentum, i.e., 
/(r,k) = /(k). To present our results, it is convenient 
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FIG. 3: (Color online). The temperature dependence of the 
pressure (a), entropy (b), and heat capacity (c). The dimen- 
sionless dipolar interaction strength are Dh ~ (dash-dotted 
lines), 0.2 (dashed lines), and 0.3 (solid lines). 



to introduce a dimensionless dipolar interaction strength 

_ ncd 

where n — N/V is the number density of the system and 
eO, = h'^{k°p f/{2m) with k°p = {Gn'^ny/^. We note that 
kp and are, respectively, the Fermi wavevector and 
Fermi energy of an ideal Fermi gas. 

We plot the phase-space distribution function /(k) 
in Fig. [T] (a) and (b) for various temperatures and 
Dh = 0.3. At low temperature, the momentum distri- 
bution is clearly stretched along z-axis. However, as one 
increases the temperature, the momentum distribution 
becomes less anisotropic. Such behavior can be most 
easily visualized by calculating the aspect ratio of the 
cloud in momentum space a = ^ {k"^) / {kD- Figure [T] 
(c) shows the temperature dependence of a correspond- 
ing to various dipolar interaction strengths. For small 
T, a approaches to the value at zero temperature; as one 
increases the temperature, a goes to unit asymptotically. 

Due to the spatial homogeneity of the system, the di- 
rect dipolar interaction energy vanishes. Therefore, the 
total energy only contains the contributions from kinetic 
energy and Fock exchange interaction. In Fig. [2] (a) and 
(b), we plot, respectively, -Ekin and i?cxc as functions 
of the temperature. The kinetic energy corresponding 
to stronger dipolar interaction is larger than that with 
weaker dipolar interaction, since the stronger dipolar in- 
teraction results in larger momentum space deformation. 
As to the exchange interaction energy, we see that i?oxc 
vanishes as a approaches 1 at high temperature limit. We 



also plot the temperature dependence of the chemical po- 
tential in Fig. [5] (c) . Clearly, /j, decreases as one increases 
T, in analogy to the ideal Fermi gas. In addition, the 
chemical potential is also a decreasing function of Dh- 
This can be most easily understood at zero t emp erature 
limit, where the chemical potential becomes [201 
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in weak interaction regime [Dh ^ 1). Equation ([6]) in- 
dicates that, for given density n, the chemical potential 
decreases as the dipolar interaction strength grows. 

To calculate other thermodynamic quantities, we work 
with the thermodynamic potential 



^l{T,V,^l) ^ E -TS - nN, 
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where 
S 
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{/(k)ln/(k) + [l-/(k)]ln[l-/(k)]}. 



is the entropy of the system. The pressure directly relates 
to the thermodynamics potential as PV = —fl{T,V,fj,) 
and the heat capacity can also be evaluated using the 
definition Cy = '^{w)vn' ^ present the 

temperature dependence of P, S, and Cy corresponding 
to different dipolar interaction strengths. As a compar- 
ison, we also plot the results for an ideal Fermi gas. In 
general, the results for dipolar Fermi gases approach to 
those corresponding to the ideal gas at high temperature 
limit. 

Figure [3] (a) also indicates that, at a given tempera- 
ture, P is a decreasing function of Dh, in agreement with 
the zero temperature result [ll[ . This can be understood 
since the overall dipolar interaction is attractive for ho- 
mogeneous gas. In addition, even though S vanishes at 
zero temperature, the entropy for given T ^ also ex- 
hibits a similar dependence on Dh as the pressure [Fig. [3] 
(b)]. Furthermore, The temperature and interaction de- 
pendence of the heat capacity, shown in Fig. [3](c), can 
be naturally deduced from the behavior of th entropy. 

To gain more insight into the entropy, we shall derive 
an analytical expression for S at low temperature limit. 
To this end, we note that when T — > 0, the derivative of 
S with respect to T at constant V and n can be approx- 
imately expressed as ^l| 
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dT 



where we have explicitly expressed / as a function of 
the temperature T and the zero temperature dispersion 
relation for the quasi-particle is [l^ 
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FIG. 4: (Color online). The phase-space distribution func- 
tion f{p,z,0,0) of a trapped gas in spherical potential for 
kBT/{N^/^hu}) = 0.2 (a), 0.5 (b), and 1 (c). The dimension- 
less dipolar interaction strength is Dt = 1. 



with P2{ ) is the second order Legendre function, 9 is the 
polar angle of k, and 



I{x) 



J3_ 3(l-x2) 



-In 



1 



l-x 



Following the same procedure as that in Ref. [21| , we find 
'dS\ V 



S{T,V,fi)~T 



dT J 
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where 
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so{k,e)=fi 



We note that Eq. ^ is not the entropy in usual sense 
as it is a function of /i. To proceed further, we change 
variable from ^ to iV by using Eq. ([B]), which yields, at 
low temperature limit and in weak interaction regime. 
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This expression clearly indicates that the entropy is de- 
creasing function of Dh, in qualitative agreement with 
our numerical results. 



B. Trapped case 

Now we turn to study the trapped case, for which we 
shall adopt a set of dimensionless units based on the har- 
monic osciUator length a = yJh/{m'Uj): N^/^a for length, 
^1/6^-1 .^yave vector, and N^/^fWO for energy. Un- 
der these choices, we may define a dimensionless dipolar 
interaction strength 



Ari/6 



Cd 



In addition, the normalization condition for phase-space 
distribution function becomes {2tt)~^ J drdkf{r,'k) = 1. 

In Fig. m we present the typical results for phase- 
space distribution function /(p, z, 0, 0) for a gas trapped 
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FIG. 5: (Color online). The temperature dependence of the 
deformation parameters a (a) and /3 (b) for A = 0.1 (solid 
lines), 1 (dashed lines), and 10 (dash-dotted lines). The di- 
mensionless dipolar interaction strength are Dt = 1. 
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FIG. 6: (Color online). The temperature dependence of the 
critical dipolar interaction strengths for various trap geome- 
tries. 



in a spherical potential. Similar to the momentum space 
distribution for a homogeneous gas (Fig. [1]), the distri- 
bution function in real space becomes elongated along 
z-axis at low temperature. However, the deformation 
becomes smaller and smaller as the temperature is in- 
creased. To characterize the real space deformation, 
we define the deformation parameter in real space as 
/3 = X^^\/ {x^) I (z^), where the inclusion of the factor 
is to eliminate the deformation caused by trapping 
potential such that P — 1 for a noninteracting gas. Fig- 
ure [5] (a) and (b) plot, respectively, the temperature de- 
pendence of the deformation parameters in momentum 
and real spaces. As expected, stronger dipolar interac- 
tion induces larger deformations in momentum and real 
spaces. However, both a and /3 approach unit at high 
temperature limit. 

Finally, we explore the stability of the trapped gases 
at finite temperature. Similar to the zero temperature 
case [1, [l^l , due to the partially attractive feature of the 
dipolar interaction, there also exists a critical dipolar in- 
teraction strength for a given temperature T, beyond 
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which the system becomes unstable. In Fig. [6l we plot 
the temperature dependence of the critical dipolar inter- 
action strength for various trap geometries. It can be 
seen that, independent of the trap aspect ratio A, 
is an increasing function of T. This can be understood 
as higher temperature corresponds to smaller momentum 
and real space deformations, resulting in smaller dipolar 
interaction. In addition, higher temperature also yields 
higher kinetic energy which stabilizes the system. 

IV. CONCLUSIONS 

To conclude, based on the semi-classical theory, we 
have studied the thermodynamics of a dipolar Fermi gas 



by numerically finding the phase-space distribution func- 
tion. We show that the deformations in both momentum 
and real space becomes smaller and smaller as one in- 
creases the temperature. For homogeneous gases, we also 
calculate pressure, entropy, and heat capacity. In partic- 
ular, at low temperature limit and in weak interaction 
regime, we obtain an analytic expression for the entropy, 
which agrees qualitatively with our numerical result. The 
stability diagram for a trapped gas at finite temperature 
is also obtained. 
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